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Fundamental to much work in small- a; QCD is a fcr-factorization formula. Normal expectations 
in theoretical physics are that when such a result is used, citations should be given to where the 
formula is justified. We demonstrate by examining the chains of citations back from current work 
that violations of this expectation are widespread, to the extent that following the citation chains, 
we do not find a proof or other justification of the formula. This shows a substantial deficit in the 
reproducibility of a phenomenologically important area of research. Since the published formulae 
differ in normalization, we test them by making a derivation in a simple model that obeys the 
assumptions that are stated in the literature to be the basis of fey-factorization in the small- x 
regime. We find that we disagree with two of the standard normalizations. 



I. INTRODUCTION 

Recently, the participants in a Les Houches meeting 
published a set of recommendations for the presentation 
by experimental groups of the results of searches for new 
physics at the LHC [lj . The recommendations are based 
on the following guiding principles: 

(El) What has been observed should be clear to a non- 
collaboration colleague. 

(E2) How it has been observed should be clear to a non- 
collaboration colleague. 

(E3) An interested non-collaboration colleague should 
be able to use and (re)-interpret results without the 
need to take up the time of collaboration insiders. 

(The labeling is ours.) 

In this paper we will argue that principles of a similar 
kind should be applied to theoretical papers, and that 
it is a substantial impediment to progress in high-energy 
physics that these principles arc not met by a large num- 
ber of published papers. We will support our argument 
with examples from the literature on the Regge, or small- 
x, limit of QCD. This area is highly relevant for the LHC, 
thereby giving a natural connection with the recommen- 
dations for experimental results. 

The processes of theoretical physics are rather different 
to those of experimental physics, so a direct transcription 
of the principles for presenting experimental results is 
inappropriate. Instead we propose the following: 

(Tl) What has been calculated or derived should be 
clear to any colleague. This has a consequence that 
for the concepts involved, either precise definitions 
should be provided or correct references to where 
correct definitions can be found. 
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(T2) It should be clear to any colleague whether the ob- 
tained results are at the level of a mathematical 
proof, a conjecture, a guess etc. 

(T3) Any interested colleague, with appropriate exper- 
tise and knowledge of the underlying theory, should 
be able to understand, use, and, most importantly, 
reproduce the results without the need to obtain 
extra explanation from the author (s). 

Now it is fundamental to science that results should be 
reproducible. In theoretical physics this means that an 
interested reader can, for example, replicate the deriva- 
tions and calculations in a scientific article. It is im- 
portant to be able to do this independently of the au- 
thor^) of the original article. The derivations in a spe- 
cific paper are not always self-contained, and references 
to where proper derivations can be found are then es- 
sential. An exception is when adequate derivations are 
genuinely part of the expected education or experience 
of the intended audience of an article. So a corollary to 
(T3) is what should be standard practice: 

(T4) When previous results are relied on, appropriate 
references must be made to where appropriate jus- 
tifications are actually to be found. This need not 
apply if the results are a standard part of the edu- 
cation of scientists working in the general area (e.g., 
theory of elementary particles). 

A general application of the principles (T1)-(T4) is nec- 
essary for the progress of our field. We consider them 
to represent the traditional standards of good theoretical 
physics. 

For reproducibility of results, there is an important dif- 
ference compared with experimental physics. Replicating 
an experimental result can often mean an amount of work 
comparable to that for the initial result, although, as time 
goes by, initially highly difficult techniques can become 
routine. (E.g., the measurement of jets in high-energy 
collisions.) In theoretical physics, the initial formulation 
and derivation of particular results may require a high 
degree of unusual creativity and insight. But the verifi- 
cation of the results should be much more routine. 
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It is not uncommon that a theoretical argument in- 
volves unstated assumptions or tacit knowledge that is 
not explicitly codified. In such cases extensive consul- 
tation with the authors of an article might be needed 
before the results of the article can be reproduced. If 
widespread, such a situation is highly undesirable and is 
in the long run inimical to good science. 

Some of the tacit knowledge may not be easy to for- 
mulate explicitly, and has to be acquired by appropriate 
guided exposition. We consider classic examples to be 
the concepts of mass and force in Newtonian mechanics, 
which resist a precise definition in terms of pre-existing 
concepts. But for the health of a science, this kind of 
tacit knowledge should be strongly limited in extent. 

To counterbalance the negative consequences of our 
finding that some results are not easily reproducible, it 
is important to remember that subjects such as those 
we discuss are fundamentally difficult, both conceptually 
and technically. It is perhaps inevitable that progress 
involves intuitive ideas acquired through long experience 
working in an area. Furthermore, it is genuinely difficult 
to convert such intuitive ideas to a fully consistent and 
self-contained mathematical structure. 

One motivation for the present paper came from an 
examination of the literature on the Regge limit and the 
related small- a; limit of QCD. These limits are widely 
applied to phenomenological studies and predictions at 
HERA, RHIC and the LHC. We noted that the concept 
of transverse-momentum-dependent (TMD), or fac- 
torization appears widely, and it is desirable to unify or 
relate the various definitions of TMD distributions, in- 
cluding those in the recent book by one of us [2, Chs. 13 
& 14]. To do this systematically requires that we know 
sufficiently accurately what the definitions are and how 
results for cross sections are obtained. But in trying to 
do this, we discovered that in many cases, comprehensi- 
ble, accurate, and valid definitions and, proofs or other 
justification, are hard or impossible to find, certainly if 
one just follows the citations given. These problems have 
some resemblance to the issues found on the experimen- 
tal side in Ref. Q, and has led us to formulate (T1)-(T4) 
above, as reasonably self-evident principles of good prac- 
tice in theoretical science. 

The related topic of the QCD formulation of the Regge 
theory is particularly interesting for our purposes. Prior 
to the discovery of QCD in 1972 to 1974, Regge theory 
was a primary topic of work in the strong interaction. 
Therefore there exists an older generation of physicists 
who are/ were knowledgeable in this subject, in particu- 
lar in the fundamentals. But with the advent of QCD, 
work naturally strongly turned away from Regge theory 
in favor of perturbative QCD. (An INSPIRE search find 
title Regge and date 1972 gave us 214 results, but 
find title Regge and date 1983 gave 18.) 

Even so, the phenomena addressed by Regge theory 
have not disappeared, and small- a; physics makes heavy 
use of it (interestingly, the search find title Regge 
and date 2010 gave 59 hits, showing an increase of in- 



terest in the old ideas). Therefore for any young physi- 
cist entering this field, there is a lot of knowledge that 
must be learned from papers that deal with rather old 
concepts whose foundation may not be particularly well- 
known today, even though the concepts are widely used. 
We are particularly concerned with foundational issues, 
such as: Why are the concepts and results in the sub- 
ject what they are? How are they defined? How do we 
know they are valid? It is then necessary to go back to 
original papers, or even simpler, textbooks on the sub- 
ject (when available). Even with textbooks and review 
articles, principles (T1)-(T4) should still apply. Q 

The issue at hand is not just that one should make 
it easier for newcomers to properly learn the field, but 
that it is also necessary to be able to verify claims of 
theoretical results made in the literature. 

Rather than give an exhaustive list of the cases where 
we have found substantial difficulty verifying results 
stated in the literature, we will take one particularly 
important case, fey-factorization as applied to the pro- 
duction of hadrons in hadron-hadron collisions. We will 
attempt to determine from some important papers on 
the subject how this particular form of factorization is 
justified. 

We emphasize that it is not necessary that there be an 
actual rigorous proof of a formula that is used. There are 
many situations in theoretical physics where it is useful to 
propose or conjecture a property or formula on the basis 
of some intuitive idea, for example, or from some natu- 
ral extension of existing ideas to new situations. Even in 
these situations it is important to know what is the jus- 
tification for a formula, and why a particular formula is 
used rather than one of the infinitely many other possible 
formulas. Therefore we prefer to use the word "justifi- 
cation" rather than "proof". Without knowing the sta- 
tus of the justification or proof and without being able 
to evaluate the arguments, it is difficult to evaluate the 
significance of subsequent work where problems are en- 
countered. 



II. AN EXAMINATION OF THE LITERATURE: 
WHERE IS THE ORIGIN OF 
fc T -FACTORIZATION? 

In this section we present explicit examples from the 
literature where we have identified problems with the ap- 
plication of the principles (T1)-(T4), in the case of fcy- 
factorization in single inclusive particle production in the 
small- cc regime. The reasons for choosing this particular 
subject are: (i) fc^-factorization has wide phenomenolog- 
ical applications, (ii) It plays a very fundamental role in 
the formulation of small-x QCD. (iii) The concept has 
been around a rather long time and therefore it is a rea- 
sonable expectation that the subject is well established 
and that our principles are obeyed, so that we can check 
the results. Any problems in formulating the theoretical 
methods can be treated as an important source of system- 
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atic error in the comparison of theory and experiment. 
Explicitly noticing such problems can be an important 
motivation for topics for further research that are impor- 
tant for the success of a field. 

The /^-factorization formula that we discuss, Eq. ([1]) 
below, is intended to be valid for single inclusive jet pro- 
duction in hadron-hadron collisions. It is widely used 
in phenomenological applications to study the particle 
multiplicity observed at hadron colliders. (For some ex- 
amples see [IHIIJ and references therein. A comparison 
of some phenomenological predictions to LHC data for 
the particle multiplicities in both proton-proton and lead- 
lead collisions was presented by the ALICE collaboration 
.13].) 

The /cT-factorization formula used in this area is (see, 
e.g., Ref. @, Eq. (1)]): 



da 



2a, 



cPpxdy CfPt 

x d 2 k A , T f A (xA,k A , T ) f B {x B ,PT ~ k A ,r)- (1) 



Here f A and fs are TMD densities of gluons in their 
parent hadrons, and the two gluons combine to give an 
outgoing gluon of transverse momentum px which gives 
rise to an observed jet in the final state. In the formula, 
C F = (N 2 - 1)/2N C with N c = 3 for QCD, y is the 
rapidity of the final-state jet, and x A: b = (pt / \fs)e ±y . 
The incoming hadrons A and B can be protons or nuclei. 

Questions that now naturally arise are: Where does 
this formula originate from? Where can a proper deriva- 
tion be found, and under what conditions and to what 
accuracy is the derivation valid? What are the explicit 
definitions of the unintegrated distributions f A ,B, and do 
these definitions overcome the subtleties that are found 
in constructing definitions of TMD distributions in QCD 
in the non-small- a; regime [3, Chs. 13 & 14]? In an ideal 
world, we could say that in order for the requirements 
(T1)-(T4) to be fulfilled, it is absolutely necessary that 
these questions be answered, and that a person who reads 
a paper which makes use of this formula can, if needed, 
go back to the original source and himself/herself repro- 
duce and verify the derivations. But we must recognize 
that in the real world some of these issues are very deep 
and difficult, and that therefore complete answers to the 
questions do not (yet) all exist. Nevertheless, in this 
subject, we should expect some kind of derivation, with 
the accompanying possibility of an outsider being able 
to identify, for example, possible gaps in the logic where 
further work is needed. 

However, as we will explain below, we tried to find any 
kind of a derivation of the formula by following citations 
given for it, but were unable to do so. Our findings can 
be visualized in Fig. [TJ which shows the chain of refer- 
ences that one needs to follow to arrive at the nearest 
possible source(s), starting from a selection of recent pa- 
pers. Coming to those sources we find that the formula is 
never derived but essentially asserted. Moreover, the ba- 
sic concepts involved are never defined in a clear enough 



way to make it understandable what exactly it is that 
is being done. We therefore find it impossible that (UJ 
can be satisfactorily re-derived from sources referenced 
in the literature, contrary to what should be the case if 
principles (T1)-(T4) hold. 

A clear symptom that these are not merely abstract 
difficulties but are problems with practical impact is that 
the overall normalization factor differs dramatically be- 
tween the references. See, for example, Eq. (40) in [l4| 
and Eq. (4.3) in [l5| — and notice that this difference in 
normalization does not appear to be commented on, let 
alone explained. The difference in normalization factors 
demonstrates that at least one of the presented factor- 
ization formulas is definitively wrong. (In the two ap- 
pendices of the present paper, we will show that in fact 
the normalizations of both formulas appear to be wrong.) 
There are a number of difficult physical and mathemati- 
cal issues that need to be addressed if one is to provide a 
fully satisfactory proof of a factorization formula. These 
issues go far beyond a mere normalization factor. But 
the existence of problems with the normalization factor 
is a diagnostic: it provides a clear and easily verifiable 
symptom that something has gone wrong. A minimum 
criterion for a satisfactory derivation is that it should be 
explicit enough to allow us to debug how the normaliza- 
tion factor arises. 

At the top of our chart of references, Fig. [TJ we have 
chosen some of the recent phenomenological applications 
[§4l"iT] that make use of ([TJ . We also include some earlier 
highly cited phenomenological applications [!, @ . There 
exist a very great number of papers which make use of 
(UJ, so we include here only a representative few. As is 
indicated in the top part of Fig. [TJ a central source that 
is given for (TJJ is the highly cited Ref. 03 . We thus ask 



whether we then can find a derivation of (QJ in 14j. 

That paper performs a calculation in a quasi-classical 
approximation of particle production in DIS using the 
dipole formalism (see the reference for the exact calcu- 
lations that define this "quasi-classical" approximation). 
There actually is an implicit assumption of a factorized 
structure from the very start in this formalism (see Eqs. 
(1) and (7) in the reference). For our purposes it is im- 
portant to notice what the exact statement is regarding 
(UJ), which can be found as Eq. (40) in |l4j. (An unimpor- 
tant difference is that in [1J], the f's in (UJ are instead 
written as //fey.) Prior to this equation, an equation for 
the production of gluons in DIS is derived, Eq. (39) in 
14j. The exact statement just prior to stating (UJ in the 
form of Eq. (40) in [3] reads 

The form of the cross section in Eq. (39) sug- 
gests that in a certain gauge or in some gauge 
invariant way it could be written in a fac- 
torized form involving two unintegrated gluon 
distributions merged by an effective Lipatov 
vertex. 

There is no derivation of Eq. (TJJ. Rather, this equation is 
stated as being the "usual form of the factorized inclusive 
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Albacete et al [10] 
1011.5161 



Levin et al [8] 
1007.2430 



Rezaeian [11] 
1111.2312 



Albacete et al [9] 
1001.1378 




Gribov et al [15] 
Phys Rep 100 
1983 



Gribov et al [22] 
Phys. Lett. B 100 
1981 



Gribov et al [23] 
Phys. Lett. B 121 

1983 



Gribov et al [25] 
Phys. Lett. B 101 
1981 



Ryskin [24] 
Yad Fiz 32 
1980 



FIG. 1. Some examples of the chain of references for the fcT-factorization formula {TJ. The boxes with crosses on top indicate 
references quoted for |TJ where the equation is not even present in the paper. The boxes surrounded by a thick red line at 
the end of a chain are where a final proof or adequate justification of the formula would be expected to be found, yet, as 
documented in the text, this is not the case. Thus the chain of references does not lead to any source from which the validity 
of the formula could ultimately be verified. Only the arXiv numbers are shown for those papers that possess one, while for 
those who do not have arXiv numbers we show the journal name and volume for the publication. 



cross section", with the source being our Refs. [1 
This is indicated by the arrows away from Ref. |l4lj m 
Fig. [TJ Furthermore, the paper then goes on to say that 
the results of the calculations actually appear to be "in 
disagreement with the factorization hypothesis" . (Our 
emphasis.) That is, Eq. (UJ is actually false in the situ- 
ation considered in Ref. |14j |. Moreover, the calculations 
in the paper are for gluon production in DIS, whereas the 
papers that cite Ref. [l4| as a source for the factorization 
formula employ it in hadron-hadron collisions. Thus Ref. 
[lij seems to be a particularly bad reference to cite for the 
validity of fcT-factorization, es pec ially in nucleus-nucleus 
scattering, as is the case in [q-lll|. 

To be fair, in a later paper the authors of [lH do 
propose a solution to the problem of the inconsistency 
of their calculations with fcT-factorization. The solution 
is that a different definition of the gluon distribution is 
needed — see the section of [3] that is titled "A tale of 
two gluon distribution functions" . But this simply rein- 
forces our statement that [14j is an entirely inappropriate 
reference to cite for fcT-factorization and its validity, at 
least not without considerable further explanation. 

We next examine the sources (l6l - [20j that Ref. 14 1 



itself cites for the factorization formula Eq. (UJ : 

• In Refs. [ljl H3| there actually exists no formula at 
all resembling (J]), and those papers are not con- 
cerned with proving factorization. We indicate this 
fact in Fig. [TJ by placing a cross on the references 

to[H,[i3. 



• In Ref. [18|, there is one equation, (40), that has 
a somewhat similar structure to our (UJ. But in 
detail it is quite different, for example as regards 
the overall normalization and the number of pow- 
ers of the strong coupling g, and the meaning of 
the functions inside the transverse-momentum in- 
tegral. Thus this reference cannot be used to sup- 
port a statement that flj is the "usual form of the 
factorized inclusive cross section" . 

• Reference [l9j is concerned with a possible interpre- 
tation of the classical formula within perturbation 
theory using the DGLAP splitting kernels, but Eq. 
UJ is never present in the reference and we there- 
fore again indicate this fact by placing a cross on 
the reference in Fig. [TJ 
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• Finally, Ref. [2(| does indeed contain Eq. (03) , which 
is referred to as the "GLR formula" and is given by 
Eq. (65) in the reference. The origin of the formula 
is given as a review article by Gribov, Levin and 
Ryskin [15j (GLR), as indicated in Fig. 00 

Thus to trace the origin of the fey-factorization formula 
(03) we have come from the recent applications down to 



the GLR review [15J, via Ref. [14|. 

So now let us turn to Ref. [Hj]. There the fey- 
factorization formula appears as Eq. (4.3). But the nu- 
merical coefficient is N c /(2tt) 6 , which is extremely differ- 
ent from the 2/Cf in our (pj) and in Eq. (40) in Q. In 
fact, as discussed in Sect. 4.4.1 of [2l|, there appear many 
different formulas for the coefficient of (Q3). The choice 
in (fTJ) corresponds to the choice in [l4{ (and the majority 
of papers use this choice, as explained in section 4.4.1 of 
21])- 

Indeed the normalization issue is referred to in the pa- 
per [2(| that is citationally nearest to the GLR paper, 
in a section entitled "Comparison with GLR formula" . 
But the combination of (65) and (68) in [20j agrees with 
our (Q3) except for a previously mentioned redefinition of 
gluon densities by a factor of fc T . No mention is made 
of the dramatically different normalization factor in the 
GLR paper, even though the GLR paper is quoted as the 
source for the fey-factorization formula. (Note that the 
formula in [2p| . like our (fTJ), is for da/cPpxdy, whereas the 
formula in GLR is for the apparently different quantity 
Eda/cPp; however, both these quantities are equal after 
applying the appropriate change of kinematic variable.) 

In addition to the issue of normalization there are sev- 
eral other relatively minor differences between the ver- 
sions of formula (00) among different papers. In [f| for ex- 
ample, the k± integral has p± as the upper limit, while 
other papers do not indicate such a limit. The use of 
such an upper limit is rather non-trivial because it does 
not treat the hadrons in a completely symmetric manner. 
Moreover, the i nteg ration measure is sometimes written 
d 2 kj_ (e.g. SSDJI) and sometimes dk\ (e.g. dUU) so 
it is not clear what the correct prescription is. Finally, 
the coupling a s is treated differently from paper to paper. 
In [lj| it is put inside the integral, with a fey-dependent 
scale, while other papers write it as in (fTJ), with a fixed 
a s outside the integral. For someone well experienced 
with this area, these differences are likely to be unimpor- 
tant. Thus the difference in the arguments of a s could 
signal the difference between a strict leading-logarithm 
approximation and an idea about the likely result of an 
improved approximation. 

But for an outsider to the topic, these differences can 
be quite bothersome, because it is not clear how they 
arise, or how significant they are, particularly when there 
is no comment on the differences. This is especially true 
in the absence of any derivations where one can locate 
the source of the differences. 

At this point, an outsider should feel entitled to be able 
to examine a detailed step-by-step derivation, so as to be 
able to pinpoint where the error (s) is/are. That there 



is an error somewhere is absolutely demonstrated by the 
difference in normalization. As we have explained, there 
are no relevant derivations in the chronologically later pa- 
pers whose citation chains led us back to Ref. [Hj]. But 
in (l5j . we find no detailed derivation either. Just above 
that paper's (4.3) there is a statement of what graphs are 
involved (with only a very brief explanation) and then a 
statement of what is supposed to be the resulting fac- 
torization formula. However, one is referred back to two 
previous papers [22|,[23[ by the same authors, as indicated 
in Fig. 03 

But in these papers, we merely find that there is an as- 
sertion that certain graphs dominate, an assertion with a 
stated reason that certain other graphs cancel, and then 
a statement of the resulting formula. There is no detailed 
derivation, e.g., a derivation that is sufficiently detailed 
to make it manifest, for example, how the many factors 
of 2ir that are ubiquitous in loop integrals organize them- 
selves to give one or other of the very different normal- 
ization factors that are found in different version of the 
fey-factorization formula. 

It should also be mentioned that the assertions made 
in Refs. [HI, [23[ are based on the axial gauge where it is 
known that the unphysical singularities introduced in the 
gauge propagators do cause problems with divergences 
in the most natural simple definitions of the TMD par- 
ton densities, and additionally the contour deformations 
needed to prove factorization are blocked by these ex- 
tra poles. See, e.g., 0, Chs. 13 & 14] for details of this. 
We find no indication in the cited references that such 
technical (yet crucial) issues are addressed. 

In [23l | , we read that the equation (with the normaliza- 
tion given in GLR, not the normalization that is currently 
accepted) "can easily be obtained in the QCD leading 
logarithmic a ppr oximation" , and the relevant references 
given are (22|.|23|. 

So we are now lead back to the papers [22|, H3] • These 
two references are indicated by the red boxes in Fig. 0Q 
and they thus constitute the two sources where one finally 
might expect to find a proof of (Q3) , but presumably with 
a different normalization factor. 

We first examine [22], in which we have been told in 
23] that a derivation is to be found. Again we find no 
derivation. There is again just a statement of the di- 
agrams involved and a statement of the resulting fac- 
torization formula. There is no derivation that can be 
debugged to check the normalization; there is not even 
any justification of why the particular graphs dominate. 
Neither is any statement given of the approximations in- 
volved. No prior references are given. Similarly in [24j 
we again find a short statement of the diagrams involved 
and of the factorization formula, but no detailed deriva- 
tion or justification. 

In [22[ it is also mentioned, below that paper's Eq. 
(1), that the authors managed to calculate the relevant 
diagrams for the high-py production of hadrons in Ref. 
[25| . We therefore also include this reference in Fig. 03 
In Ref. [25j | . the fey-factorization formula is presented in 



6 



the very last equation of the paper (un- numbered). It 
is, however, again simply stated without any derivation. 
Nor is any further reference or argument provided, and 
it is not at all clear what the definitions of the parton 
distributions are. It should also be noted that in [25| . 
an axial gauge is used for studying the leading graphs 
of high-pT processes. Thus one would here again expect 
the problems mentioned above that stem from unphysical 
gauge singularities. We find again no indication in the 
references that these problems are dealt with. 

In Fig. [TJ we have enclosed the boxes [HI, [24], (25[ in 
thick red lines. This denotes that they are at the end 
of the chains of citations, and that the derivations we 
should expect to find in them do not exist. 

In the appendices of the present paper, we present 
a derivation of a fc^-factorization formula starting from 
what appear to be the same assumptions and approxima- 
tions used by GLR. We get a quite different formula, as 
regards the normalization. Our proof gives what we in- 
tend to be enough detail that an outsider can verify the 
result. The assumptions and approximations are valid 
in a model that uses lowest-order graphs, so that our 
derivation is adequate to test whether a claimed result is 
actually true. 

A summary of this section is that we find that regard- 
ing the important and widely used factorization formula 
(UJ), the principles (T1)-(T4) are violated very badly. A 
confusing chain of citations is present in the literature as 
shown in Fig. [TJ By following the relevant references we 
are unable to find any paper where a proper derivation 
is provided. Moreover, the explicit definitions of the ba- 
sic entities involved are not always clear, and differ from 
paper to paper (see also section 4 of [21j for a more ex- 
plicit discussion) . It is therefore extremely important for 
progress in the field that all of recommendations (Tl)- 
(T4) are adapted and maintained universally. 



III. SUMMARY 

Our principles (T1)-(T4) embody sensible standards 
and conventional norms for the presentation of work in 
theoretical physics. Suppose a paper uses, but does 
not derive, a particular formula of the status of the kr- 
factorization formula (UJ as a key part of its work. Then 
it is a normal and reasonable expectation that a refer- 
ence should be given to where the formula is derived or 
otherwise justified. This we found to be badly violated, 
as summarized in Fig. [TJ 

For example, a main reference fl4| frequently given for 
(TJJ not merely failed to derive it but actually questioned 
its validity in the relevant context. Moreover, of the ref- 
erences that [3] gives for the formula, most failed to 
even contain the formula. Backtracking through a chain 
of references where the formula does appear we eventu- 
ally met a dead end: we could not actually find a useful 
derivation. 

This is a quite disturbing situation. The formula in 



question is a key part of an important area of QCD phe- 
nomenology, and a superficial reading of the literature 
indicates that the formula appears to be accepted with- 
out question as the standard one. But closer examination 
fails to turn up a justification, at least not by following 
the normal and uncontroversial procedure of examining 
the references cited for the formula. 

Of course, it may be that by a much wider search 
through the literature one can find an adequate justi- 
fication of (JJ). But the point of a citation is to avoid 
a broad literature search: The reader just needs to go 
to the cited article(s) and find both the formula and its 
justification. Any substantial deviations from this expec- 
tation need explicit explanation. 

A failure to find appropriate justification of a main- 
stream formula immediately raises the question as to 
whether the formula is correct. Recent work 26 28] finds 
violations of standard factorization in related situations, 
so ignoring the problems is both dangerous and can lead 
to failure to find correct consequences of QCD. 

For a single paper whose citations are inadequate, it 
may be most appropriate to contact the author(s) to 
get corrected information. But in the present case, the 
problems are much more widespread, and need a more 
widespread concerted remedy. We also emphasize that 
we have identified other similar problems in the jus- 
tification of related topics but for the sake of brevity 
we document here only the situation for one case: fcy- 
factorization for gluon production in hadron-hadron col- 
lisions. 

We conjecture that one reason for the inadequate ci- 
tations is that the /cy-factorization formula is regarded 
by its users as uncontroversial. Thus providing cita- 
tions for it is a formality rather than a serious scientific 
exercise. (29j We conjecture that a second reason is that 
there is substantial pressure, particularly on young sci- 
entists, to produce new results, and that it is regarded as 
counterproductive to actually check the validity of earlier 
results that are relied on to obtain the new results. 

In addition, a disincentive to giving a full derivation 
is its length, if one makes explicit the details needed 
by a newcomer. For example, the derivation in the ap- 
pendix of the present paper adds about 50% to the pa- 
per's length. 

Our findings show that such tendencies need to be re- 
duced for the health of our field; they greatly and unnec- 
essarily impede the ability of other physicists to repro- 
duce and verify the work. Insiders may have good reason 
to trust the formula, but workers in even closely related 
areas of QCD cannot check the results, at least not at all 
easily. The responsibility is not only on authors to justify 
their results adequately, but also on journal referees to 
check that results quoted in a new paper have actually 
been justified. Note that the problems we are finding are 
not at the level of a deep critique of a difficult proof; we 
find there is such an absence of derivation that there is 
not even a proof to critique. 



7 



ACKNOWLEDGMENTS 

This work was supported in part by the U.S. Depart- 
ment of Energy under grant number DE-FG02-90ER- 
40577. We thank Markus Diehl, Ted Rogers, and Mark 
Strikman for very useful comments on a draft of this pa- 
per; the current version includes responses to some of 
their questions. 



Appendix A: Derivation of the normalization factor 
of the fey-factorization formula 



Since we have been unable to hnd a derivation of the 
fey-factorization formula ([I]) by following references for 
it, and since there is disagreement on the normalization 
factor, we provide here a simple derivation that gives the 
normalization. 

The derivation starts from certain assumptions con- 
cerning the graphs and momentum regions that domi- 
nate. Since the assumptions are demonstrably incorrect 
at sufficiently high order in perturbation theory, as we 
will explain, the derivation cannot be completely correct 
as regards full QCD. Nevertheless, the derivation does 
apply to a simple Feynman graph model, and therefore 
unambiguously determines the appropriate value for the 
normalization, which differs from both of those in the 
standard references, [H, Eq. (4.3)] and [14, Eq. (40)]. 
Moreover, the assumptions we make are compatible with 
the stated starting point given in [l5|, [2^, l23l j , so our 
derivation can be used to test the correctness of the fac- 
torization formula stated in those references. 

Our derivation is a simplified version of a derivation 
that one of us gave in Ref. |21j . 



1. Definition of process and gauge 

We consider the cross section da/cPprdy for inclusive 
production of a gluon in a high-energy collision of unpo- 
larized hadrons, where the gluon has transverse momen- 
tum pt and rapidity y. The kinematic region of interest 
is where A <C \pt\ *C ^fs and e VB <C e v <C e VA , where 
y/s is the center-of-mass energy, and yA and ys are the 
rapidities of the beam hadrons. We define Q = \pr\, and 
let the momenta of the incoming hadrons be Pa and Pb ■ 

The aim is to obtain a factorized form for the cross 
section that is valid up to corrections that are a power of 
the small parameters pr/ ' \fs, e~( y ~ VB \ and e~( VA ~ v \ 

We use light-front coordinates relative to the two 
beams, with the momentum of the detected gluon being 



,P ,Pt) = 



Qe v Qe 



V2 V2 



-,Pt 



(Al) 



use the usual symbol A for another purpose. The gauge- 
fixing vector has the same rapidity as the detected gluon: 

noc (e v ,e- v ,Q T ), (A2) 
Then the numerator of the gluon propagator has the form 



Wn" + n^k u V l k»? 



k ■ 



(k- 



(A3) 



Following [30j,[31|, we use a time- like axial gauge n-<j) = 
0. Here the symbol <fi denotes the gluon field, since we 



The point of this gauge is that for gluons of high positive 
rapidity it approximates the light-like axial gauge <f> + = 
while for gluons of high negative rapidity it approximates 
the light-like axial gauge (\f~ — 0. This happens since 
k ■ n ~ k + n~ when the rapidity of k is large and positive, 
while k ■ n ~ k~n + when the rapidity of k is large and 
negative. 

When various complications are ignored, each of these 
different light-like gauges gives an appropriately simple 
formulation of parton-model physics inside the corre- 
sponding beam hadron — see 0, Sec. 7.4] and references 
therein. 



2. Assumptions 

We make the following assumptions: 

1. The cross section is dominated by (cut) graphs of 
the form of Fig. [5J where one gluon out of each 
hadron combines to form the observed gluon. Lines 
associated with the incoming hadrons are in the 
subgraphs A and B. 

2. We assume that the rapidities of lines in A are sub- 
stantially greater than the rapidity y of the ob- 
served gluon, i.e., the rapidities in A are at least 
y + A, with e A ^> 1. This and the next assump- 
tion encode a key kinematic property used in BFKL 
physics. 

3. Similarly the rapidities in B are less than y — A. 

4. Transverse momenta in A are of order Q or smaller, 
and transverse momenta in B are of order Q or 
smaller. 

5. We ignore the fact that a massless gluon is not a 
physical particle in real QCD, as opposed to low- 
order perturbation theory. 

6. We ignore all other subtleties, e.g., that certain 
other graphical structures may be important, and 
can even break factorization |26l428| . 

The assumptions are valid for the lowest-order graph 
in a perturbative model of high-energy quark-quark scat- 
tering, Fig. [3] In this graph, it is readily checked that 
the kinematic assumptions for the A and B subgraphs 
follow from the kinematic region chosen for the detected 
gluon. It can also be checked that other graphs of this 
order are power-suppressed in the gauge we have chosen, 
given that polarization vectors e for the final-state gluon 
obey e • p = e ■ n = 0. 
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FIG. 2. Graphs assumed to dominate for Eq. (0. 



/ 



Pr 



Pa 
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FIG. 3. Lowest order example of Fig. [2] Here n, a, etc, 
are Lorentz and color indices for the gluons attaching to the 
horizontal gluon line. 



Cross section 



At high-energy, the cross section from Fig. [5] is 



da 



1 f d A k A d 4 k B 



d 2 pTdy 2s 
1 



(2tt)4 



5 {i) {k A 



p)x 



16tt 3 



j^lfl 1 W -qvv' ,(}[}' jj i 



(A4) 



Here, A MM ■ aa is the amplitude for the subgraph A in- 
cluding the propagators for its two external gluons, as a 
function of the Lorentz and color indices of these gluons. 
It includes integrations over internal loop momenta and 
the momenta of the final-state particles of A. Similarly, 



B vu',pp' is for subgraph B. The factor 

1 



16tt 3 



(A5) 



arises from the Lorentz-invariant differential for the de- 
tected gluon: 



d 3 p dyd 2 p T 



{2irf2E p 16tt 3 ' 



(A6) 



given that the cross section is differential in px and y. 
Finally, H^i , V u';aa' ,pp' is the factor for the central gluon: 



P lifi' ,vv' ;aa' ,/3/3' ^ ^ 9 fafi~ffa 



'/3' 7 - 



x [(k A - k B ) ■ eg^ v + (k B + p) M e„ + (-p - k A ) v e l _ l ] x 

x[(k A - k B ) ■ eg^'u- + (k B +p)^e v > + {-p - k A )^e^} . 

(A7) 

The sum is over two independent orthonormal polariza- 
tion vectors for the gluon. In a frame in which px is along 
the x-axis: 



/ e v e -y \ 



(AS) 



we can choose 

e a = (0, 0, 0, 1) and e 2 = (e y , -e~ v , 0, 0) /V2 (A9) 

as the independent polarization vectors. 

To keep track of the dependence of calculations on 
the gauge group, we will assume that the gauge group 
is SU(iV c ). 

4. Leading-power approximations 

We now make approximations suitable for the kine- 
matic region that we have assumed. 

First we need the size of the components of the mo- 
menta k A and k B of the gluons coming out of A and 
B: 

k A = (9^. + o(Qe v - A ), o(Qe-y- A ), k A>T ) , (A10) 
Qe~y 



k B = 0{Qe 



v/2 



+ 0(Qe- y ~ A ), 



(All) 



The components kT^ and kg are small relative to the cor- 
responding components of p, since the components of k A 
and k B cannot exceed the typical values for momenta in 
A and B; the sizes of momenta in A and B are governed 
by our assumed separation of rapidity between A, B and 
the gluon p. Then momentum conservation at the three- 
gluon vertex ensures that k\ and kg equal p + and p~ , 
up to small corrections. 
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Next consider the line for one of the external gluon 
lines of A. It connects a subgraph H characterized by 
rapidity y to a subgraph A with rapidity y + A, and we 
write the relevant part of the graph as 

a x N^(k A )h». (A12) 

Here is the numerator of the gluon propagator, while 
a and h correspond to the rest of the A subgraph and to 
the H subgraph. We treat /i M and a x as being obtained by 
a boost from a situation with zero rapidity, and therefore 
we obtain relative sizes for the components: 

a = {a+,a-,a T ) ~ const, x (e v+A , e~ y - A , 1), (A13) 
h ~ const, x (e v , e~ v , 1). (A14) 

We use the sizes given in (|A10[) to determine that the 
sum over components in (|A12[) is dominated by A = + 
and the transverse components in fj,: 

- a+N-' l K l . (A15) 

A corresponding calculation is then applied to each of the 
gluons connecting the central gluon to the B subgraph. 

Finally we observe that k A and k^ differ from Qe y j\[2 
and Qe~ v /y/2 by an amount that is small compared with 
the corresponding components of momenta in A and B. 
Thus in A and B we can replace k A by Qe y / v2, and kg 
by Qe~y/V2. 

We then find that, given our assumptions, the cross 
section is given by 

da ir f 2 r dkT^ r dkg 
d^prdy = 2s J J (2^7 (2^* 

x ,aa' gjf ,p/3' jjii' ,jf ;aa' ,p/3' (A16) 

where the integrals and factors of 2ir are arranged to 
correspond to those in the definition of a gluon density, 
as we will see. 



where 4>^' a denotes the gluon field operator. This defini- 
tion is equivalent to Eq. (2.3) of [30( when the light-cone 
gauge <j) + — is used, which we will now show to be 
appropriate in the approximation we are using. 

If gluon momenta all have high rapidity, then k ■ n = 
k + rT [\ + |qr ~ k + n~ in the gluon propagator, 

where the neglect of k~n + j (k + n~) relative to unity fol- 
lows because the rapidity of k is much higher than that 
of n. Note that for the gluon momentum k A entering 
the hard scattering, the size of the plus component is 
appropriate for the hard scattering, i.e., k A ~ Qe v , but 
the minus component is appropriate for subgraph A, i.e., 
^a ~ Qe~ v+A , so that k A has effectively rapidity of 
about y + A. 

Hence in the gluon density, time-like axial gauge is 
equivalent to <f> + = gauge to leading power. The above 
definition then gives 

P(x,k T ) = k+ J —A u > aa (k,P). (A18) 

The normalization of the above definition is that 
P(x, kr) would be a number density of gluons in x and 
kx, if we ignored the complications present in higher or- 
der corrections in QCD. The corresponding quantity in 
the fcT-factorization formula ([T]) is actually a momentum 
density. Therefore, in the above formulas, we follow Refs. 
[3CJ, |35[, and use P(x, kr) to denote the number density, 
rather than the more common notation f(x, kx)- 

Of course, as explained in 0, Chs. 13 & 14] the light- 
cone-gauge definition [3B| has problems, because of rapid- 
ity divergences and other issues, and even the definition 
with a non-light-like axial gauge has trouble. So the def- 
initions have to be modified to be suitable for full QCD. 
But within the approximation used here, the simplest 
definition is adequate. 

6. Factorization 



5. Gluon density 

If there were no complications, the natural method to 
define a number density of a parton is by the hadron 
expectation value of the number operator in the sense of 
light-front quantization. The basic ideas were worked out 
by Bouchiat, Fayet, and Meyer [32[ and by Soper (33|, 
although these papers do not contain actual formulas 
for transverse-momentum-dependent parton densities in 
terms of hadron matrix elements of field operators. [34] In 
light-front quantization and light-front gauge, this gives 
the gluon density as 



P(x,k T ) 



xP^ 



duu d wt e 



(2tt) 3 

x (P|^ a (0,u;-,^r) 



We next recall that in the small- a; limit and in the 
approximation we use here, only the term (IA15|) in the 
lines for the external gluons of A is important. So for 
each gluon line there is a factor of the vector transverse 
momentum multiplying a rotationally invariant quantity. 
Moreover with a color-singlet or a color-averaged beam, 
the dependence of A on color is proportional to 8 aa ' ■ 

Thus for the A factor in (IA16|1 . we have 



dk A 



-A 1 ' 



dk. 



\k AtT \ 2 N*-lJ (2tt 

A A °aa' 



\k 



A,T 



N? 



lfc+ 



P(x A ,k AtT ) 



(A19) 



Here the factors of k\ 



and k A arise from the axial- 



|P), (A17) 

gauge gluon propagator as it is used in (IA15|l . while 
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xa = Qe v / (^/2p\) — Qe v /y/s (in the overall center-of- 
mass frame). An exactly similar result applies to the B 
subgraph, with x B = Qe~ v / (y/2pg) — Qe~ v /y/s, Sub- 
stituting these results in (|A16j) gives 

^ , = 7T2 / d 2 k A , T P A {x A ,k A ,T)x 

d z p T ay x A x B s z {N* - ly J 

XPB(x B ,PT-k A , T f AkA ^ B " . (A20) 

\kA,T\ \ k B,T\ 

Performing the algebra in the H factor, and making ap- 
proximations valid to leading power in e _A then gives 
the following factorization formula 

da 2w 2 a s f 2 

-T, T = ~F< T / " KA,T XaPa(xAiKA,T)X 

d z p T dy CfPtJ 

x x b Pb{x b ,Pt - k A)T ), (A21) 

with a s = g 2 /{4n) and C F = (N 2 - 1)/(2JV C ), as usual. 
Notice that each gluon density has been multiplied by a 
factor of the corresponding x value, and that this com- 
pensated a change from 1/s to \/p\ in the pre-factor. 

This factorization formula is to be compared with the 
formula stated earlier, ([lJ, which has the normaliza- 
tion given in [1J, Eq. (40)]. Evidently in that formula, 
the gluon densities are to be interpreted as momentum 
densities rather than number densities, i.e., /(ie, &t) = 
xP(x, fcr). But there is also a factor of tt 2 missing. 



7. Why the proof is insufficient 

Our derivation of the factorization formula (|A21[) relied 
on certain assumptions, which are not in general true in 
QCD. However they are true for the lowest-order graph, 
so that our derivation gives the correct normalization fac- 
tor. At a minimum the proof needs enhancement if the 
result is to be regarded as actually derivable in full QCD. 
We now explain some of the difficulties. 

In the first place the definition of the gluon density 
needs to be changed. Already, it was seen in [13, HH 
that the natural definition of a TMD quark density in 
light-cone gauge has rapidity divergences. This was cor- 
rected by the use of non-light axial gauge. But even this 
has problems [37], App. A]. A satisfactory gauge-invariant 
definition for a TMD quark density was recently given by 
one of us in 0, Chs. 13 & 14], and the definition can be 
extended in a reasonably natural way to a TMD gluon 
density. However, the treatment in [2| did not cover the 
additional issues that arise at small x. 

In general, many more kinds of graph and momentum 
region contribute. If factorization is to hold, various sum- 
mations and cancellations have to be employed. More- 
over some of the cancellations involve issues of relativistic 
causality — see [HI, [3(| — which is broken graph-by- 
graph by the k ■ n denominators in the axial-gauge gluon 



propagator — see (|A3|) . Use of the Feynman gauge in- 
creases the number of graphical structures and regions to 
be considered. 

Furthermore, it has recently been found that factoriza- 
tion is actually violated in some significant situations in 
hadron-hadron collisions [2ol-|28|. While these particular 
examples do not include the particular process discussed 
in the present paper, they do show that claimed proofs of 
factorization need careful attention to ensure that they 
actually work. 



Appendix B: Testing the normalization factor 

The argument that derives the normalization of the kr- 
factorization formula is somewhat long. Since the result 
disagrees with both of the standard formulas to be found 
in the literature, it is useful to find some elementary tests 
to verify our result. 

In our formula (|A21I) the factors of tt 2 are different 
from those in both [H Eq. (4.3)] and 114, Eq. (40)]. The 
group theory factor 1/Cp does agree with that in [Tij . 
but disagrees with the factor N c in [l5| , which is written 
there simply as N. 

We do the verification by determining how the pow- 
ers of tt and the N c dependence arise in the lowest-order 
graph Fig. [3] for the cross section and in the lowest-order 
graph for the gluon density. We will show how this im- 
mediately gives the power of tt and the N c dependence 
in the prefactor in the factorization formula (|A21[) . In 
the lowest-order graphs and in n ■ <f> = gauge, all our 
assumptions are valid, and agree with the starting point 
in GLR [l5|, so at lowest-order there is no ambiguity as 
to the correct result. 



1. Power of 7r 

In the cross section, the integral over the plus and mi- 
nus components of the loop momentum are performed 
using the mass-shell delta-functions for the spectator 
quarks. After the standard kinematic approximations 
these give simple algebraic calculations, and hence no 
factors of tt. There remains a two-dimensional trans- 
verse momentum integral which is unchanged between 
the cross section calculation and the factorization for- 
mula. From the differential (|A6|) , there is a 1 /tt 3 for each 
of the 3 final-state particles, and there is a tt 4 associated 
with the delta-function for overall momentum conserva- 
tion. This gives an overall factor of 1/tt 5 in the cross 
section. 

In the one-loop gluon density the definition itself con- 
tains an overall factor of 1/tt 4 , and there is a factor of 
tt from the 2tt8((p — k) 2 ) that puts the spectator quark 
on-shell. This gives an overall factor of 1/tt 3 . 

To match these powers of tt, the coefficient of the two 
gluon densities in the factorization formula for the cross 
section has a factor of tt. Since a s = g 2 /(47r), we need a 
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factor 7r 2 multiplying a s , which is exactly what we have 
in (|A2ip . The factorization formulas in both of Refs. 
flil Il5l] have a different power of n and are therefore 
wrong. 



Here repeated indices are summed, as usual, and the 
result was obtained by use of the formulas tr(£ Q t Q /) = 
\&aa', fapffafs'y = N c 6pp>, and <5 77 = N 2 — 1. 

Each parton density has a group-theory factor 



2. Dependence on N c 

We now ask for the dependence on N c when the gauge- 
group is SU(iV c ). We choose the quarks in Fig.|3]to be in 
the fundamental representation, and we choose the cross 
section to be averaged over the colors of the initial-state 
quarks. In a correct factorization formula, the coefficient 
is independent of the nature of the beam particles, of 
course, so these choices for the beams are purely for cal- 
culational convenience. 

Then the group theory factor for the graph for the cross 
section has the form 



1 



(Bl) 



1 N 2 - 1 

tv(t a t a )=C F = ^— . 



(B2) 



The coefficient in the factorization formula therefore has 
a factor 



N?-l 1 



AN C C 2 N 2 -l 2Cp 



(B3) 



in agreement with our factorization formula in (|A21[) and 
with the one in [l4j , but in disagreement with [l5j. 



[1] 

[2] 
[3] 

[4] 
[5] 

[6] 

[7] 



S. Kraml, B.C. Allanach, M. Mangano, H.B. Prosper, 
S. Sekmen, et al, "Searches for new physics: Les Houches 
recommendations for the presentation of LHC results," 
(2012), [arXiv:1203.2489 [hep-ph]| 

J. C. Collins, Foundations of Perturbative QCD (Cam- 
bridge University Press, Cambridge, UK, 2011). 
Mark Strikman in a personal communication remarked 
that it is an interesting question as to whether the liter- 
ature on Regge theory meets these standards sufficiently 
well. But we do not wish to answer that particular ques- 
tion here. 

Dmitri Kharzeev, Yuri V. Kovchegov, and Kir- 



sion in pA collisions," 


Phys.Rev. D68, 094013 (2003) 


arXiv:hep-ph/0307037 [hep-ph] 





Nestor Armesto, Carlos A. Salgado, and Urs Achim 
Wiedemann, "Relating high-energy lepton-hadron, 
proton-nucleus and nucleus-nucleus collisions through 
geometric scaling," |Phys.Rev.Lett. 94, 022002 (2005) | 
|arXiv:hep-ph/0407018 [hep-phj| 

Dmitri Kharzeev, Eugene Levin, and Marzia 
Nardi, "Color glass condensate at the LHC: 
Hadron multiplicities in pp, pA and AA col- 
lisions^ iNucl.Phys. A747, 609-629 (2005) | 

|arXiv:hep-ph/0408050 [hep-phj| 

Eugene Levin and Amir H. Rezaeian, "Gluon 
saturation and inclusive hadron production 
at LHC," |Phys. Rev. D82, 014022 (2010) | 

arXiv:1005.0631 [hep-ph]' 



[8] Eugene Levin and Amir H. Rezaeian, "Hadron mul- 
tiplicity in pp and AA collisions at LHC from the 
color glass condensate," |Phys. Rev. D82, 054003 (2010)] 
|arXiv:1007.2430 [hep-phj| 
[9] Javier L. Albacete and Cyrille Marquet, "Single inclusive 
hadron production at RHIC and the LHC from the color 
glass condensate," |Phys.Lett. B687, 174-179 (2010)] 
arXiv: 1001. 1378 [hepjphl 



[10] 

[11] 
[12] 

[13] 

[14] 

[15] 
[16] 

[17] 

[18] 
[19] 



Javier L. Albacete and Adrian Dumitru, "A model 
for gluon production in heavy-ion collisions at the 
LHC with rcBK unintegrated gluon densities," (2010), 
|arXiv:1011.5161 [hep-ph]] 

Amir PL" Rezaeian, "Charged particle multiplicities 
in pA interactions at the LHC from the color 
glass condensate," Phy s.Rev. D85, 014028 (2012), 
|arXiv:llll 2312 [hep-ph]] 

Prithwish Tribedy and Raju Venugopalan, "QCD satu- 
ration at the LHC: comparisons of models to p+p and 
A+A data and predictions for p+Pb collisions," (2011), 
|arXiv:1112.2445 [hep-ph]] 

B Abelev et al. (ALICE Collaboration), "Charged- 
particle multiplicity density at mid-rapidity 
in central Pb-Pb collisions at ^/ s jy at = 2.76 
TeV," Phys .Rev.Lett. 105, 252301 (2010) | 



|arXh7: 1011.3916 [nucl-exj 
Yuri V. Kovchegov and Kirill Tuchin, "Inclusive 
gluon production in DIS at high parton density," 
|Phys. Rev. D65, 074026 (2002) | |arXiv:hep-ph/0111362| 
L. V. Gribov, E. M. Levin, andM. G. Ryskin, "Semi hard 
processes in QCD," |Phys. Rept. 100, 1-150 (1983)| 
Alex Kovner, Larry D. McLerran, and Herib- 
ert Weigert, "Gluon production from non-abelian 
Weizsacker- Williams fields in nucleus-nucleus 
collisions," |Phys.Rev. D52, 6231-6237 (1995)] 

|arXiv:hep-ph/9502289 [hep-ph] | 

Alex Kovner, Larry D! McLerran, and Heribert 
Weigert, "Gluon production at high transverse momen- 



structure functions," 


Phys.Rev. D52, 3809-3814 (1995) 


arXiv:hep-ph/9505320 [hep-ph] 





Yuri V. Kovchegov and Dirk H. Rischke, "Clas- 
sical gluon radiation in ultrarelativistic nucleus- 
nucleus collisions," |Phys.Rev. C56, 1084-1094 (1997)| 
|arXiv:hep-ph/9704201 [hep-ph] | 

Xiao-feng Guo, "Gluon minijet production in nuclear col- 



12 



lisions at high-energi es," |Phys.Rev. D59, 094017 (1999) 
|arXiv:hep-ph/9812257 [hep-phj] - 

[20] M. Gyulassy and Larry D. McLerran, "Yang- 
Mills radiation in ultrarelativistic nuclear col- 
lisions/^ |Phys.Rev C56, 2219-2228 (1997) | 

|arXiv:nucl-th/9704034 [nucl-th]p 

[21] Emil Avsar, "TMD factorization and the gluon 
distribution in high energy QCD," (2012), 
|arXiv:1203.1916 [hep-ph]j 

[22] L. V. Gribov, E. M. Levin, and M. G. Ryskin, 
"High pt hadrons in the pionization region in QCD," 
|Phys. Lett. B100, 173-176 (1981)| 

[23] L. V. Gribov, E. M. Levin, and M. G. Ryskin, "Large- 
Et processes as a main source of hadrons at very high- 
energies," |Phys. Lett. B121, 65-71 (1983) 

[24] M.G. Ryskin, "Hadron productin with large transverse 
momentum in the pionization region and the vacuum sin- 
gularity in QCD," Sov. J. Nucl. Phys. 32, 133-140 (1980), 
yad. Fiz. 32, 259-274 (1980). 

[25] L.V. Gribov, E.M. Levin, and M.G. Ryskin, "Extend- 
ing the QCD leading logs up to reggeon field theory," 
|Phys.Lett. B101, 185 (1981) | 

[26] John Collins and Jian-Wei Qiu, "kr fac- 
torization is violated in production of high- 
transverse-momentum particles in hadron- 



hadron collisions," 


Phys. Rev. D75, 114014 (2007) 


arXiv:0705.2141 hep-pl 





[27] Ted C. Rogers and Piet J. Mulders, "No 
generalized TMD-factorization in hadro- 
production of high transverse momentum 



hadrons," [Phys.Rev. D81, 094006 (2010) 


arXiv:1001.2977 [hep-ph] 





[28] Jeffrey R. Forshaw, Michael H. Seymour, and Andrzej 



Siodmok, "On the breaking of collinear factorization in 
QCD," (2012), [arXTv:120676363 [hep-ph]| 
[29] In addition, the foundational work on a topic, which may 
be decades old, may give the key ideas, but not the extra 



steps to put results in the form in which they are later 
found to be most useful and convenient. It may be that 
these extra steps are known privately to insiders but that 
they do not appear in the papers that it is most natu- 
ral to cite. We give an explicit example in Sec. IA 5l This 
phenomenon can explain the problems with the citations, 
but it leaves unchanged the readers' difficulty of repro- 
ducing results. 

[30] John C. Collins and Davison E. Soper, 
"Parton d istribution and decay functions," 
|Nucl. Phys. B194, 445 (1982)| 

[31] John C. Collins and Davison E. Soper, "Back-to-back jets 
in QCD," |Nucl. Phys. B193, 381 (1981) | 

[32] C. Bouchiat, P. Fayet, and P. Meyer, "Galilean invari- 
ance in the infinite momentum frame and the parton 
model," Nucl. Phys. B34, 157-176 (1971). 

[33] D. E. Soper, "The parton model and the Bethe-Salpeter 
wave function," Phys. Rev. D15, 1141-1149 (1977). 

[34] This gives an example of the difficulty referred to in foot- 
note [B2] 

[35] D. E. Soper, "Partons and their transverse momenta in 
QCD," Phys. Rev. Lett. 43^ 1847-1851 (1979). 

[36] Note that the definition in [2j is for a TMD quark density. 
For a TMD gluon density, one makes the same adjust- 
ments as are used to get the formula for an integrated 
gluon density from the formula for an integrated quark 
density. 

[37] Alessandro Bacchetta, Daniel Boer, Markus Diehl, 
and Piet J. Mulders, "Matches and mismatches in 
the descriptions of semi-inclusive processes at low 
and high transverse mom entum," |JHEP 08, 023 (2008) 
|arXiv:0803.0227 [hep-ph]| 

[38] G. T. Bodwin, "Factorization of the Drell-Yan cross- 
section in perturbation theory," Phys. Rev. D31, 2616- 
2642 (1985), erratum: D34, 3932 (1986). 

[39] J. C. Collins, D. E. Soper, and G. Sterman, "Soft gluons 
and factorization," Nucl. Phys. B308, 833-856 (1988). 



